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1 Introduction 

In this paper we develop extensions of the conical techniques exposed in the 
book and we develop new algorithms as well. We make an effort to keep 
this paper enough self-contained, however, the book ^S, (hereafter often referred 
to as "the book" for the sake of brevity) can be useful for further details. Such 
book is an outgrow of two papers published on the journal Optimization, namely 
ID and m. 

To begin with the conical approach has his own nomenclature. Terms like 
duality have a different meaning than those of standard approaches (see e.g. 

Other terms like strict tangency are peculiar to the conical approach. The 
term bounded refers to a LP problem (maximization of a linear functional on 
a polyhedron) that is feasible and has maximum. The specification internal 
for a primal algorithm has also a meaning related to the conical context, and, 
as explained below, refers to the particular approach followed here to reach 
optimality, as opposed to the primal algorithm presented in the book. 

In the book we stressed the particular interest of conical conditions that are 
at one time purely conical ( that is, involve only cones) and purely pointed (that 
is, the involved cones are pointed). Dual conical feasibility and optimality con- 
ditions are of this kind. Although they have lead in the book to an enumerative 
algorithm, the dual conical methodology is not only useful for its peculiar prop- 
erties (e.g. it gives a closed form expression of the maximum, solves a whole 
class of problems and is particularly efficient in handling parameter variations), 
but it has provided an useful tool for applications in fields like control and game 
theory. The related literature is by now very extended and relevant, and has 
provided many interesting new results (see e.g [5] and the bibliography therein). 

As to the primal conical conditions given in the book, the first is conical but 
not pointed. The second involves a pointed cone and an affine space, and has 
been the starting point for the development an evolutive primal external conical 
algorithm, enjoying exact finite convergence. 
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We shall take here this whole business to a more advanced stage. First be- 
cause we shall give a further primal condition, that is purely conical and purely 
pointed. Secondly because, intertwining this new primal approach with the gen- 
erators machinery used in the book to develop the dual conical algorithm, we 
provide a further algorithm, which although based on generators is evolutive. 
Thirdly we give a closed form expression of the maximum in the primal setting 
too. More than that we give explicit closed forms of the solution both on the 
range space of the coefficient matrix and in terms of unknowns of the LP prob- 
lem. Finally we stress that the present result complete the theory presented in 
the book from yet another point of view. In fact the second primal optimality 
condition is a tangency condition (of an afhne space to the non-negative orthant) 
and the previous primal algorithm reaches tangency landing the afhne onto the 
cone. Thus we knew that one could in principle try to reach tangency from the 
other side: that is, starting with the two sets that meet each other and then 
taking the affine to touch the only boundary of the cone (thereby achieving tan- 
gency). However, developing such an algorithm (emerging submarine instead of 
landing hydroplane metaphor) has defeated us for a while. It is finally presented 
here. 

As to notations we stick to those adopted in the book with only one variation 
in the interest of simplicity. The variation is that when we shall get to the pa- 
rameterized feasibility formulation of optimality we shall denote the augmented 
coefficient matrix and bound vector in the same way as the non-augmented case, 
leaving to the context the specification of which is which. Recall also that we use 
always the same symbol P for the non-negative orthant of any euclidean space. 
Again the space to which such symbol is referred is specified by the context. 

2 Mathematical base 

We will assume thorough the hypothesis of strict tangency of the feasibility and 
linear programming problems we study. Whereas this is no restriction of gen- 
erality as proved in the above book, it has a great geometrical importance and 
it yields a much more neat and elegant path to the development of algorithms. 
Note that actually such an assumption is made in the primal conical LP al- 
gorithm introduced in and reported in the book, within the theorem that 
states exact finite convergence. Here however, we simplify matters avoiding to 
postpone the use of such hypothesis as much as possible in an effort to achieve 
maximum generality. 

Let us start considering a linear feasibility problem. That is, the problem of 
determining whether or not a polyhedron is non-void and in the positive case 
finding at least a point in the polyhedron. Better yet would be finding more 
points, i.e., having the possibility of exploring the polyhedron, even, ideally, 
finding all of its points. A polyhedron is the intersection of a finite set of 
semispaces. Thus a set Q of the form: 

g = {x:Gx <v} 
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where G is an n x m real matrix, x (the unknown vector) is in i?™ and v 
(the bound vector) is in i?". The polyhedron is a cone if and only if v = O.We 
shall often denote "^(G) (the range of G) by F, for brevity. Recall also that the 
vector V — Gx is called the slack vector (a feasible slack vector \i Q ^ (f) and 

X e g). 

We now make a simple but important remark, independent of the hypothesis 
of strict tangency. If we decompose v as v = vj? + where v^? is the orthogonal 
projection of v on and v is the orthogonal projection of v on _F , we can 
write for the inequality defining the polyhedron: 

Gx <WF+v 

Because vector inequalities are translation invariant, assuming Gz — "vp, the 
latter is true if and only if: 

Gx — = G{x — z) < V 

Thus in passing from one to the other condition, feasibility is invariant, the 
slack vector too, and the solution is varied by a constant vector z. Note that 
the polyhedron is a cone (when v^- = 0) or a translated cone (when vp ^ 0). 
if and only if w = 0. In what follows we assume that this is not the case i.e.: 
V =/= 0. Otherwise we always have the trivial solution a; or, respectively, 
x — z, and our investigation would become pointless. Also notice that under 
this assumption the slack vector cannot be zero because v Ww € 
R^,v -Gw T^O. 

Remark 1 Notice that if v E P, then Q ^ (f), because x = Q is obviously a 
solution. Similarly ifvEP, then Q ^ (p, because x = z is obviously a solution. 

To say that the associated problem is strictly tangent |3] means that the 
subspace 7^(G) is strictly tangent to the non-negative orthant P, that is, Tl{G) 
meets P in the only origin: 

n{G) n p = {0} 

We now recall the second primal conical feasibility (necessary and sufficient) 
condition |2|: 

g ^ ^ (v + TZ{G) nP^(j) 

This is a conical condition, but there is an affine set appearing in it, namely 
v + 7^(G). 

Recall also that, if feasibility prevails, to derive a feasible slack vector in the 
range space of G we can find any vector y in the intersection: 

y e (v + 7^(G) nP 

Consequently, if we want a solution in the domain space, we can solve in x 
the system: 
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Gx = V — J/ 



Obviously if we solve this equation for all y in (v + 7?.(G) fl P we get all and 
nothing but the points of the polyhedron Q. 

3 New primal conical feasibility conditions 

We shall now start elaborating further on this second primal feasibility condi- 
tion. One of the dividends will be the introduction of a new purely conical and 
pointed primal condition. As a first remark note that (as it is immediate to 
prove): 

V + 7^(G) = v + n{G) 

and therefore we can write: 

Next we consider a subspacc that is extended with respect to Tl{G) , namely 
(with £(.) we denote linear extensions and with Co(.) conical extensions): 

Fe = £(v) +F = C{v) +F = C{v + n{G)) 
Finally we need the cone: 

Ce = Co{v) + F = Co{v + F) 

Note that clearly: 

Fe = Ce + (-Ce) 

We are now in a position that allows to state the following new primal conical 
feasibility condition, which requires that a pointed cone do not degenerate into 
the trivial (singleton of the origin) cone: 



Theorem 2 The following primal conical feasibility conditions hold: 

g^cp^Ce n P 7^ {0} 

or equivalently: 
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Proof, li Q ^ cf), by the above version of the second primal feasibihty 
condition, {v + 7?.(G) D P ^ (j) and this intersection is made up of non-zero 
vectors. Let y be one of those, so that y = v + z with z in TZ{G). Therefore 
there arc non-zero points in Ce H P and the condition is necessary. 

Conversely take y 7^ in Ce n P. It will have the expression y = fiv + z with 
z in 'R.{G) for some real /3 > (if /? were zero y would be zero too by virtue of 
strict tangency). Thus the non-zero vector w = l/l3{f3v + z) satisfies: 

w = 1/PiPv + z) = v+^zeCe 

But evidently it is also true that w € {v + 'R-{G)) fl P and this completes the 
proof. ■ 

The procedure explained in the proof to pass from the vector y in the cone 
Ce n P to the vector w moving along the ray generated by y, will be called 
calibration. Notice that such procedure is viable numerically in a very simple 
way. In fact, if we denote by Ppj. the orthogonal projection onto F-^, it is 
obvious that: 

p = {PF±y)i/vi;yi s.t. Vij^O 

Another important fact that follows immediately from the above proof is 
that if we let, in the feasibility case Pc = (i" -|- 7?.(G) n P (more on this set later) 
we can state the following: 

Corollary 3 The following relation holds whenever feasibility prevails 

Co{P,) = Ce n P 

Next notice that if Q = (f) so that Ce fl P = {0} it may either be the case 
that: 

-Ce n P = {0} 
in which case it is clear that Pe is strictly tangent to P; or that 

-Ce n P 7^ {0} 

In this respect, the next natural question is to find out whether it can be the 

case that — Cg fl P {0} and Cg fl P 7^ {0} at the same time. We shall show in 
the proof of the next result that this cannot be the case, because the two cones 
Ce and — Ce live in opposite semispaces and the non- negative orthant, with 
the exception of the origin, is entirely contained in the interior of one of them. 
And this means that we can establish the following further feasibility condition, 
which is more handy in that it substitutes the subspace Pe to the cone Ce- 
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Theorem 4 If D P ^ {0} then the sign of (3 is constantMy ^ in F^nP. 
Moreover: 



Fe n Pj^ {0} and P>0 



Proof. Suppose F^ O P ^ {0} and consider any ^ in F^ H P. Because F 
is strictly tangent to P, we know from theorem 6.2.1 in 3 that F is contained in 
an hyperplane H , which is strictly tangent to P and there is a vector n, normal 
to this hyperplane, which is internal to P. Then for some z 7 ^ : 



Assume to fix the ideas that 7 > so that (n, u) > 0, y G Ce and therefore 
y S.Ce n P. In this case by the previous condition feasibility prevails. Clearly 
by the same argument and because n is internal to P: 



Thus both the cones Ce and P are contained in the semispace {x : (n, x) > 
0}, delimited by the hyperplane H, and, in addition, P is in the interior of the 
semispace, with the only exception of the origin, because n is in the interior of 
P. 

Moreover:: 



and this latter implies that —Ce is in the opposite semispace (i.e. the semis- 
pace {x : (n, x) < 0}) with respect to P, so that —Ce Ci P = {0}. Therefore the 
sign of /3 in Fe n P = (Ce + (-Ce))n P is constant. 

Assuming instead that 7 < and reasoning along the same lines (n, v) < 0, 
y £ — Ce and therefore y G.— Ce fl P. Moreover in this case: 



and therefore Ce H P = {0} and feasibility does not prevail. In this case 
— Ce n P ^ {0} and the sign of /3 is constant in — Ce n P. Thus the proof is 
finished. ■ 

Corollary 5 There are three mutually exclusive and exhaustive cases possible 
for the feasibility problem: 

a) Fe is strictly tangent to P. That is Fe D P = {0}; In this case the 
problem is unfeasible. 

b) Fe r\ P ^ {0} and /3 < 0. In this case too the problem is unfeasible. 

c) Fe Ci P ^ {0} and f3 > 0. In this case the problem is feasible 



< {n, y) = {n, {^v + z)) = 7(71, v) 



(n,Ce) C [0,+oo) 
{n,P) C [0,+oo) 



(n, -Ce) C (-00,0; 
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As in all algorithms exposed in once it is determined that the problem 
is feasible and an ?/ G Ce n P is found, a solution in the domain space can be 
found solving the equation Gx — v — w, where w is obtained from y by the 
calibration procedure. 

4 Feasibility algorithm 

We next introduce a feasibility algorithm, based on the results obtained so far. 
To this purpose we have to devise a method to find (in the feasible case) one or 
more feasible solutions in the range space y £ Ce H P . Or, which is more easily 
viable, find an y G nP with (3 > Q and then calibrate it. The method should 
fail if and only if unfeasibility prevails so to completely solve the problem. One 
possible and useful way to go is to exploit the generators machinery for pointed 
polyhedral cones (and we got purely conical and pointed conditions primarily 
to that effect) in this primal conical setting. The advantages of this approach 
will be more and more evident in the sequel. Because y must be in a pointed 
polyhedral cone given by the intersection of a subspace with the non-negative 
orthant, we can find one or more of the extreme rays of this cone or, equivalently, 
of its generators. A good part of the book is devoted to the development of 
a sophisticated machinery to solve this problem, with basic algorithms and two 
levels of more advanced algorithms as well as an implementation and numerical 
results. The fundamental results on which the generators techniques are based 
are given in Chapter 9. 

To carry out this program, we only need to find an operator whose null space 
is Pe- This is not difficult. In fact: 

Pe - C{v) + P 

Where the sum is a direct orthogonal sum. We introduce the notations V = 
L(v) and Pv and Pp to denote, respectively, the orthogonal projections onto 
V and P. Then the orthogonal projection onto C{v) + P is given by Py +Pf 
so that / — Py — Pp is the orthogonal projection onto {V + P)^ Therefore: 

Pe = C{v) + P = N{I -Pv- Pf) 

Thus we can state the following feasibility algorithm. In the feasible case 
one can compute either a single or more solutions in the range space. This is 
not a detail it is a major point on which we shall expand in the next section. 

Algorithm 6 ; New primal conical feasibility algorithm 

1- Preliminary step: verify whether v or v are in P. If either of those is 
the case trivial solution are immediately available as explained in Section 2. In 
such a case STOP. 

2- Compute Py and Pp and hence I — Py — Pf ■ 

3- Use the techniques of '31 to compute a first generator of J\f{I — Py — Pf)^ 
P. If none is found the problem is unfeasible. Hence STOP. Otherwise go to 
step 4- 
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4- Calibrate the generator. (Recall that (3 = cannot happen) If f3 < the 
problem is unfeasible. In such case STOP. If (3 > the problem is feasible. 
Proceed to the next steps. 

5- Facultative step. Can be executed or skipped. Co on to calculate more or 
all the generators and calibrate each of them. 

6- For each calibrated generator gc a solution can be obtained solving Gx = 
V- 9c. 

Remark 7 In certain applications, it may be convenient ( of course in the fea- 
sible case) to find a point in relative interior of the pointed polyhedral cone 
FeCl P = Af{I — Py " Pp) n P. To this effect one can compute all the calibrated 
generators and then the sum of them (Theorem 6.1.2 in [3J) is in the relative 
interior of the cone. 

The significance of the set of all the calibrated generator goes way beyond 
this simple remark and is of crucial importance in the present context. This 
issue is fully explored in the next section. 

5 The contact polytope 

Suppose that the problem is feasible. Then, as recalled above, the polyhedron 
= (v + 7^(G) n P = {v + TZ{G) n P is non-void. In the book, leaning on 
a result by Ben Israel it is shown that this set is actually a polytope, which 
we called the contact polytope. The properties of the contact polytope play 
an important role in the book, but it was not determined explicitly. Here we 
complete the theory with a detailed description. In fact the set {gd : i = 1, .., k} 
of calibrated generators of the pointed polyhedral cone Ce n P is exactly the 
set of extreme points of Pc. This is stated in the next theorem: 

Theorem 8 The set of calibrated generators of the cone Ce C\ P is equal to the 
set of extreme points of the contact polytope Pc = (v + TZ{G) H P. Thus: 

Pc - C({.gc,}) 

Proof. First we prove that any ray of the cone G^ C) P can intersect P^ only 
in a single point. For suppose that two points of a ray, say zi and Z2, be in 
the intersection, so that we may assume, without restriction of generality, that 
Z2 = OLZ\ with a > 1. Then Z2 — zi 7^ would be both in TZi^G) and in P. But 
this contradicts the hypothesis that Tl{G) is strictly tangent to P. Therefore 
Z2 = zi as we wanted to prove. 

We know from the previous section that all the points of {gd : i = I, ..,k} 
are in Pc and Pc is contained in Ce H P. Let now {zj : j = 1, ..,p} be the set of 
extreme points of Pc. 

We start showing that each calibrated generator is an extreme point of Pc. 
In fact suppose that is not so for a certain g^. Then we can write (in what 
follow we drop indexes in sums to keep the notation simple): 
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with at least two non-zero coefficients. Isolating one term, we can write without 
restriction of generality (possibly the sum reduces to a single term): 



where a\ + P = 1 and w £ Pc- Because z-i is an extreme point, zi ^ w, and 
by what wc proved at the beginning, these two vectors cannot be proportional. 
It follows that this equality contradicts that lies on an extreme ray and 
therefore we have reached the conclusion that all the points Qd are actually 
extreme; points of P,.. 

Finally assume that there is some extreme points that are not in the set of 
calibrated generators. Let one of these be Zj. Then because Zj.€ Ce n P: 



where ai > V i and ^ > If just one a, is non-zero a contradiction 
is immediate, because if cti = 1, then Zj is a calibrated generator and if 7^ 1 
there are two proportional vectors in Pc- Excluding this case, either X^aj = 1, 
and then we get the contradiction that an extreme point is a non-zero convex 
combination of a set of other extreme points or J2 7^ 1- If that is so, dividing 
both sides of the above expression of zj by ^ a, we get again the contradiction 
that there are two points in a ray that belong to the polytope Pc- The proof is 
thereby completed. ■ 

The determination of the contact polytope and, in particular, the identifi- 
cation of the extreme rays of Cg n P with the rays generated by the extreme 
points of Pc is a very important addition to our theory. We can score on this 
result immediately. 

We start with a remark on solutions. In fact Pc can be also viewed as the 
set of all feasible slack vectors, in the sense that it contains all and nothing but 
the slack vectors y for which all solution of Gx + 2/ = v is a feasible solution. 
Consequently the set of vectors x obtained in this way is the set of all feasible 
solutions. We can summarize this in the following- 
Corollary 9 if we let the above feasibility algorithm run to compute all the 
calibrated generators, then all and nothing but the solution of the problem are 
given by {x : Gx = v - y, y e C{{gci})}. 

Remark 10 It should be stressed at this point that the present theory yields in 
a primal conical setting an explicit expression for both the polytope of feasible 
solution in the range space (slack vectors) and the polyhedron of the feasible 
solutions of the problem 

Another important consequence of this result is that we are now in the 
position of introducing an internal primal conical LP algorithms. 



9ci = aizi + ^ 



a, 



Zn = aizi + an) X] '^^^^'> = "^^^ + ^'^ 
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Before getting into this (in the next section), let us recaU briefly the well- 
known parameterized feasibility formulation on which we based all LP methods 
(see the book for more details). 

Consider the problem of maximizing a linear functional on a polyhedron 
(Problem LP) 

max fx 
subject to : Gx < v 

We can rewrite this as: 

max h 
subject to : Gx < v(/i) 

where the augmented coefficient matrix G is obtained adding to G a row with 
the entries of — / and the augmented bound vector v(/i) is obtained adding to 
V a last entry equal to —h. 

In what follows we shall soon use for this problem the same notations of the 
feasibility problem in order to simplify our presentation. Thus the reader is ad- 
vised that it is the context to determine whether we refer to plain or augmented 
coefficient matrix and to plain or augmented bound vector. The same rule we 
apply to all the other mathematical entities related to the problem. Thus, for 
example we still denote by Pc the contact polytope corresponding to the poly- 
hedron defined by Gx < v(/i). Also the non-negative orthant is still denoted by 
P in the augmented euclidean space. However, whenever beneficial to clarity, 
we explicitly denote dependence of the appropriate items on h (writing, e.g., 
Pcih)). ^ 

In this parameterized feasibility setting, to solve the problem, we have to 
find an ho such that if h > ho then the polyhedron {x : Gx < v(/i)} is void and, 
if ft, < ho, the same polyhedron is non-void. Such an ho, if it exists at all, is the 
optimum value of the functional. In terms of the second primal conical feasibility 
condition we have to find ho that verifies the following tangency condition: 

(v(/i) + 7^(G)) nP = (f>iih>ho 

(v(ft) + n{G)) r\P = Pc^4>\ih<ho 

The primal conical LP algorithm in the book (under conditions that insure 
feasibility and boundedness - see next section) started from a large value of h 
, so to insure that Pc be void and diminished such an h until the tangency of 
the affine space to the non-negative orthant is reached, so that such value of h 
is just the optimum value hooi the functional . This is the external approach. 

Here; we take the dual (how many meaning of this word!) view: we start 
from low values of h, so that the contact polytope is non-void, and increment h, 
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until the contact polytopc is squeezed into the maximal face of the non-negative 
orthant, that corresponds to a zero last component. At that point tangency 
prevails and, therefore, the corresponding value of h is equal to the optimum 
ho- Actually we give an algorithm of this sort in two versions. The second one 
will realize a further particularly important advance within our methodology. 
To each of them we devote the next two sections. 

6 Primal conical internal algorithm: first ver- 
sion. 

As in the book, and without restriction of generality, we assume strict tangency 

Feasibility can be ascertained as illustrated in the previous section. We recall 
from the book that, in view of Theorem 6.3.1, under strict tangency, feasibility 
implies boundedness. Thus we can now work on optimality with all three hy- 
pothesis (strict tangency feasibility and boundedness) in force. Keep in mind 
that we refer to the augmented problem now, although we do not change nota- 
tions. 

In view of the results of the foregoing section, in the feasible case, to know 
the calibrated generators of Ce n P is the same as knowing the extreme points of 
the contact polytope Pc- As we shall see momentarily the knowledge of certain 
extreme points of Pc. is equivalent to the knowledge of the maximum value of 
the functional . It cannot be overemphasized the importance of the fact that, 
by the following theorem, we obtain a closed form expression for the maximum 
in a primal setting, just as a closed expression for the maximum was already 
given in the dual setting (Theorem 5.3.1 in the book). 

Theorem 11 Suppose to choose an h such that the augmented problem is 
feasible(h < ho) and let {gd : i = 1, .., k} be the set of extreme points of Pc{h) 
(or, what is the same, calibrated generators of .Ce{h) (1 P) then, letting hm = 
max{fi(ci„+i : i = 1, ..,&}.• 

ho = h + hm 

Proof. Passing from an /i to /i -|- A/i corresponds to adding — A/i to the 

last component of all points of v + TZ{G). Thus all points in the contact 
polytope with the last component equal to hm are still in the new contact 
polytope Pc{h + hm) and fall in the maximal face of P, M„+i = {y : y € P and 
yn+i = 0}. Suppose now that in Pc{h+h,-,i) there is a point z with Zn+i = p > 0. 
Then passing from h+ hm to h (that is, incrementing h by —hm) this point z 
translates to a point with all the same components but the last, which is equal 
to p + hm > hm and, in addition, it clearly is in Pc{h). However the last 
component 7 of all the points in Pc{h) satisfies hm = inaxjgcj^^^ : i = 1, .., k} 
> 7 > min{(/ci„+i : « = 1, .., k} by fact that a polytope is the convex extension 
of the set of its extreme points. Therefore we have found a contradiction and it 
follows that Pc{h + hm) C M„+i, or, in other words, v -|- Tl{G) is tangent to P. 
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By the second primal conical optimality condition, ho ^ h + h,n is the optimum 
value of the functional and solves the problem in the range space of G. ■ 

In the proof of the theorem the procedure to find solutions is already built 
in. We record such procedure in the following Corollary (by e,; we denote the 
vector that has all zero components but the ith, which is equal to 1): 

Corollary 12 (i) Let gc be an extreme point of Pc{h) such that 5c„+i = hm- 
Then an optimum slack vector is given by Ho — gc" hmCn+i, Ho is an extreme 
point of Pc{ho) and all solutions of the equation Gx — v — yo (that necessarily 
exist) are optimum solutions. 

Let {gci} be the set of extreme points of Pc{ho) ■ Then the set of all optimum 
solutions is obtained solving the equation Gx = v — j/ where y G C{{gci})- 

Remark 13 With this Corollary we have completed the picture of the conical 
approach giving also the explicit closed form for the sets of optimal solutions 
both in the range space and in the domain space. 

The proof of the Corollary requires but trivial new verifications ans can be 
safely omitted. 

We can now structure a PL algorithm.. First set an h small enough (it can 
be arbitrarily small) to ensure that h < ho {a, remark on this is given right after 
the statement).. Then apply the following 

Algorithm 14 (Primal conical internal algorithm) 

Stepl Find all the calibrated generators with positive last component. 

Step2. Set ho = h + h„i , where hm is the maximum of last components of 
the found calibrated generators 

StepS Consider any calibrated generator g^ such that 5cj„+i — h„i. Then an 
optimum slack vector is given by yo — gd ~ hmCn+i and an optimum solution 
Xo is given by any solution of the equation Gx = w ~ yo. 

The proof of the algorithm is given in the theory so far developed. The 
above Corollary also illustrate how to find the set of all the solutions, if needed. 

Remark 15 Note that the internal and external algorithm (Described in Ch 11 
of the book) complete each other. Ln fact if no calibrated generator with positive 
last component is found we are either at the optimum i.e., h — ho, or h > ho. 
In any case we can revert to the external algorithm ]^ and find the solution. 

An interesting aspect of this algorithm is that it is not purely enumerative 
because we do not look for all the generators. We may easily reformulate it 
requiring that in STEP 1 all the generators be found. In this case the remark 
is changed accordingly. If we find generators, but none has a positive last 
component then /i = /iq. If no generator is found h > ho and we can revert to 
the external algorithm. 

We can pursue that feature further and introduce an evolutive version of 
the algorithm. Although the evolutive character could be exhibited in abstract 
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terms, it becomes more evident if we take to the fore .the techniques illustrated 
in the book to find the generators of a polyhedral pointed cone, which is the 
intersection of a subspace and the non-negative orthant. We conjugate in this 
way the generators technique used to deploy the dual conical methods (which 
were essentially enumerative) with the present new primal approach and ful- 
fill the quest for evolutiveness mentioned in the book within the dual conical 
framework. 

7 Computation of generators and the evolutive 
version of the algorithm 

We assume the same hypotheses and, in particular, feasibility and h < ho Our 
first purpose is to show that it is possible to apply to the present problem 
the machinery developed in the book for the computation of generators. The 
peculiarities of our method will then allow us to derive an evolutive algorithm. 

Recall that the orthogonal projection Pp^ of the space onto Fg = C{v) + F = 
V + F, is given by: 

Pf, =Pv + Pf 

and the orthogonal projection of the space onto F^ is given by: 

Pf^^ = I - {Pv + Pf) 
It follows that we can express F^ as: 

Fe = //{Ppl^ ) - AA(/ - [Pv +Pf)) 

At this point wc can apply all the machinery developed in the book to find 
the generators of 

Fe n P = Mil - (Pv + Pf)) n P 

Notice that in this formula only Py depends on h. 

However, with respect to the case of the dual conical method, there are 
numerous simplifications. First because we want those generators that have a 
non-zero last component. In this respect we can state the following: 

Theorem 16 Under the present hypotheses, deleting the last column of the 
matrix I — {Pc{v} + Pf), we obtain a matrix with the same rank as the original 
matrix. 

Proof. In view of theorem 15.1.1 of the book, if it were not so, we would not 
get any generator with non-zero last component, and hence in view of Theorem 
above a contradiction would arise. ■ 
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Corollary 17 In the procedure for search of generators of M{I —{Pv +Pf))^P 
with non-zero last component given in the book we can constantly use the last 
column as test column. 

We are now ready to introduce an evolutive conical algorithm 

Algorithm 18 (Primal conical internal evolutive algorithm) 

Put /i(0) = h. Put T(0) = I — {Pv{h) + Pp)- Repeat the following step for 
1^1,2,..: 

STEP i: If i=l perform the procedure of the hook modified fixing the last 
column as test column to find generators of .Af{T{0)) n P. If i>l resume the 
search from the sequence of basic column subsequent to the last of step i-1. If a 
generator is found, do not verify it was already found, and proceed to calibrate 
it. Let g{i) he the calibrated generator. Set h{i) = h{i — 1)+ g{i)\n+i. Compute 
T{i) setting h ~ h{i). 

Until the procedure introduced in the hook terminates. 

If the loop is exited ati = j, set hg = h{j). Optimal solutions can be obtained 
solving for x the equation G{x) = v{ho) — 9{j)- 

Theorem 19 The evolutive algorithm enjoys exact fijiite convergence. That 
is it converges in a finite number of step and if the loop is exited at i = j, 
ho = h{j). 

Proof. The proof is essentially contained in the proof of the enumerative 
version. What we do is to pass from a calibrated generator to the next with 
increasing last component. Because the number of calibrated generators is finite 
it is granted the algorithm converges in a finite number of steps. The only 
thing that remains to be proved is that the technique to find the sequence of 
calibrated generators is correct. But this too is rather obvious. In fact each 
time we increase h the current calibrated generator (as well as any calibrated 
generator with lower last component) is eliminated. Thus any new calibrated 
generator has a larger last component with respect to the former ones, and, 
consequently, we never have to verify that we find already known calibrated 
generators. Moreover we can resume the search from where it left at each step, 
because if wc started from scratch and found a calibrated generator, the same 
calibrated generator would have appeared before, by the argument used in the 
proof of the first version of the algorithm, and would have appeared with a lower 
last component. And wc know from the theory developed hitherto that this is 
a contradiction. The rest of the algorithm (computation of solutions) should by 
now obvious. ■ 

Remark 20 Notice that, as made clear by the above proof, the algorithm will 
usually get rid of some of the generators that are computed in the non-evolutive 
version. In other words evolutiveness is not just adjourning the value of the 
maximum, hut, in general, avoiding the necessity of visiting the whole set of 
extreme points of the contact polytope. 
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8 Conclusion 



As for the previous conical algorithms we deferred submission until we had 
evidence that the algorithm performs correctly numerically. During the devel- 
opment of the implementation we used the same example of [J. Of course, as 
in the previous cases, timeliness was priviledged, so to arrive to a first straight- 
forward implementation, that in the present case was written in Pascal, within 
the Delphi environment.. 

The previous experience showed that optimization of the code is a lengthy 
and painful endevour, that required the derivation of further results that are 
accounted for in the book. However, we got the divident of entire orders of 
magnitude improvings in computing time. The same process is in its inception 
for the present algorithm. An uprise of its numerical efhciency will be given as 
soon as we will feel that the level optimization of the code is satisfactory, and 
further improvements will have marginal effects only. 
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